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Abstract 

We calculate various Wilson loop averages in a pure SU (iV)-gauge theory on a 
two-dimensional sphere, in the large limit. The results can be expressed through 
the density of rows in the most probable Young tableau found in |Q] . They are valid 
in both phases (small and large areas of the sphere). All averages for self- intersecting 
loops can be reproduced from the average for a simple (non self-intersecting) loop 
by means of loop equations. 
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1 Simple non self-intersecting loop 



Two-dimensional Yang-Mills theory is quite a trivial example of a quantum field theory, 
since it has no propagating degrees of freedom. Nevertheless, the study of different gauge 
invariant quantities shows quite a rich structure and can be viewed as a useful tool for 
the verification of mathematical hypotheses about the corresponding theory in physically 
relevant dimension; say, for the search of a string theory describing multicolour QCD. 

As was demonstrated in P, Q on the example of Wilson loops, and in P] on the 
example of partition functions on the compact two-dimensional manifolds, the results 
can be represented in terms of a sum over minimal coverings ("soap films") spanned 
on a contour, or wrapping over the closed manifold. In this respect, the theory can 
have a beautiful mathematical application, as was the case with the application of three- 
dimensional Chern-Simons gauge theory. 

First we shall consider a simple loop on a sphere, separating areas Ai and A2, with 
which a holonomy U is associated. The average over all possible gauge fields inside and 
outside the loop gives the partition function of corresponding discs with fixed holonomies 
U and around the boundary, each given by the heat kernel (on the group) 



{l\e-^^\U) = Y.dRe-^'n (1) 

R 

W{A,,A2) = (^Trf/) = 1| dU^TTU{l\e-^^\U){U\e-^^\l) (2) 

where A denotes the laplacian on SU{N), with characters as eigenfunctions: 

{Axr){U) = c)^ xr{U), c\i being the quadratic Casimir of the irreducible representation 

R. 

The character expansion of W{Ai, A2) reads: 

WiA,,A2) = I E ^1^2^ / f/t/Trt/xi(t/)^^e-w^i+#^2 (3) 



and 



Ri,R2 



where / dinrUxiiU)x2iU) is the number of occurrences of R2 in fundamental (8> -Ri, 
that is: 1 when the Young diagram of R2 is obtained by adding one box to the diagram 
of Ri or otherwise. 

Now, if the rows in the first diagram have lengths (?T.i)j=i,,,Ar, the corresponding di- 
mension is 

d= n (4) 

l<i<j<Af 3 ^ 

and the Casimir 

TV 

c = Y^ni{ni + N + l- 2i) - N {n)^ (5) 

i=l 
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with (n) = J2j If the i*'^ row can be added a box, the dimension is multiphed by 

4 ^ -pr j -nj-i + Ui + l 

and the Casimir changes according to 

c,/iV = c,/iV + 2(^ + l-M) (7) 

We use the notation $j = i — y — rij + (ra) and ^{x) its large N hmit: (j){i/N) = ^i/N; 
we have to compute 

W{A, A,) = 1 E ^ E n (1 + e-^- e^^^(VA^) (8) 

where the sum over i corresponds to all possible ways of adding one box to the diagram. 
This is the expectation value 

(^En(l + ^^)e^^'^^/^^) (9) 

computed through averaging over all possible Young diagrams with the same weights as 
those which appear in the partition function of a sphere with total area Ai + A2] but in 

has been found the shape of the most probable diagrams: if x corresponds to a line to 
which a box can be added, we have {Ai + A2)(l){x) = 2 dy-^^^^p^^ and we can solve for 
= obtaining a semi-circle lawQ if A = Ai + A2 < n'^ or an expression involving 
the complete elliptic integral of the third kindQ ii A = Ai + A2 > tc'^ (see |]T[]). Indeed, 
as the same weight is given to any representation and to its conjugate (represented by 
n'i = 2{n) — riN+i-i), the unique solution for (j) is symmetric: 0(1 — x) = — 0(a;), and 
corresponds to (pseudo-)real representations. 

In a similar way, we can average over R2 and sum over all possible ways of suppressing 
one box, so that VF(v4i, A2) is also given by the average 

ij^T^Ua + ^^re-^^^^^n (10) 
We turn to the computation of 

where the $j's are distributed according to the above-mentioned law, with fluctuations 
of order 1. When j is far apart from i {\j — i\ of order A^), the logarithm is of order 1/A^, 



■^The even density is = ^j^-\/(a^ — (/'^)(0^ — &^)ni( — ^5 ^ = f ) when b<(j)<a or p=l when 

U'2 i 



-b<(j)<b, with aA = AK , a{2E - k''^K) = 1. 
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and there are N such terms: these give a finite contribution to the sum. But the same 
is true for those j's that are close to i: there are a few such terms, each of order one, 
and one has to take them into account, adding their contribution to the (principal value) 
integral which accounts for the terms of the first kind: — f f^,, . . Note that the latter 
terms were discarded in the computation of the shape of the most probable diagrams: 
looking for a maximum of ci^exp ( — ^) with its continuous formulation really means 
looking for a diagram with maximal weight with respect to long-wavelength perturbations 
{5ni — 5ni+i = but for few i's, far apart from one another), in which case the contribution 
we are going to consider vanishes; while here it does not, because we consider a short- 
wavelength perturbation 6ni = 5i^i^. So, the solution of |T| has to be understood in this 
way, and we have to consider the contributions of all the "most probable" diagrams, which 
differ from one another through short-wavelength fluctuations. 

As the integers have independent fluctuations of order 1, it seems reasonable to 
replace — $j (for j — i oi order 1) by its continuous approximation N{(f){i/N) — / N)) , 
which is: {i — j)(j)'{i/N) in the large A^ limit. Thus, the contribution to the sum of 
logarithms is 

-EE^(^-jrV(0) (12) 

If we introduce f{x) = J2k>i ^x^''{B2k/i^2k)\) {Bp = p^^ Bernoulli number), and use the 
characterization of the most probable distribution, we obtain: 

i 

Noting that xf'{x) = x/2 — 1 + x/(e^ — 1), we have f{x) = ln(^^^) — | and 

W{A,,A,) = ^^Y.^^'^l^'^^ (14) 
N ^ up 

with 5 A = A2 — A1] starting from an average over R2, and considering all possible ways of 
suppressing one box, we would obtain the same expression with A2 and Ai interchanged: 
indeed, the sign of SA does not matter as the even distribution p has only even moments. 

What about the sum over all acceptable lines? Deflning n{x) to he rii/N {x = i/N), 
we see that if |^| > 1, generically we have rii > nj+i and all lines are likely to be extended. 
However, if |^| < 1, there is a flnite fraction of lines for which rii = nj+i, and the non- 
increasing condition on the n^'s prevents a flnite fraction of lines from being extended. In 
this case, the contribution of — ^ in the range dx is (ia:!^! = c?x(l — 1). 

acceptable i 

Finally, we would expect: 



WiA, A2) = I rf0Min(- - 1, l)e'^/' ^^^^ (15) 

J P TT 



1 1 ,^JA/2 </,sin(7rp) 

However, as Boulatov flrst noticed 0, one has to be careful in expressing the product 
YljjjLi{l + ^.^^. ) in terms of the continuous function 0, when |^| < 1. In this case, he 
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groups factors in packets and obtains a continuous approximation to 62/ di which is the 
same as the preceding one, but for a factor that compensates for the Min(. . .). Indeed, 
when < 1, fluctuations of lengths of the rows are prevented by the non-increasing 
condition on the rij's: this is the very reason why a direct continuous approximation is 
not valid. We shall rather consider the heights of the columns, which are free to fluctuate, 
so that the argument given above for lines gets transposed. 

When < 1, generically we have : rii — nj+i = or 1, so columns are labelled by 
n; summing over all possible additions of one box means summing over all n's. Just as 
in the preceding case, the change in the Casimir is ^^^^ = —20; as for dimensions, the 
contribution of j's far apart from i reduces to the exponential of the same principal value 
integral, recognized as exp ( — ^^^"^^ 0). To study the contribution of j's close to i we 
consider a truncated diagram, formed with a flnite number of lines extending above and 
under let denote the height of the n*^ truncated column {vi < n < 1/2) and set 
= n — hn , (j) = ^/N . How does the dimension change when we add a box at the end 
of line that is at the bottom of column n = rii + 11 

Before the extension, Y{j^i{^i — $j) contains all integers ranging from 1 to /i„ + z/2 — Ti + 1 
but the hn — h,j + p — n + 1 , z/ = n + 1, . . . , 1/2- Similarly, the product Wjyii^j — 
evaluates to {^n-^u,-'2)\ 



v\<u<n 



After the extension these factors are respectively ( j;" + ^ 2 — g^j^^ 



n i^^-^n) n (<^'n-$.)' 

n<u<U2 ui<u<n 

SO that: 

d2 T-r / 1 \ - - 1 / N 

(we have just expressed the (change in) dimension in terms of heights of columns for 
"steep" diagrams: rii — nj+i = or 1). 

The last factor goes to 1 for large truncated diagrams, extending symmetrically around 
i (use hp — hn = ~ ^) for large A^, and note that this symmetry corresponds to 



the cancellation of C(2p + 1) in (|T2|) and to the identification of the other terms with a 
principal value integral). 

Finally, arguing that the heights of columns are integers with fluctuations of order 1, 
and can be approximated by a continuous function, we apply the same computation as 
before to get: 

W{A,,A2) = ]-Y.e'^'''''^^ (17) 
where p=^ = — % = 1— p. This is exactly: 



W{A,,A2)= Idc^e^^'^^'^^ 

J TT 

The same result was independently obtained by Boulatov 
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2 Other loops 



Before showing how the loop equations |^ allow us to compute self-intersecting Wilson 
loops, we cast ([l^) into a more convenient form for later purposes. 



Define G{z) = J . 



if X hes in the support of p and is such that p(x) < 1 (a; 



corresponds to a region where boxes can be added) we have 

G{x ± it) = —X =F mp{x) 
with A the total area; if p{x) = 1 then sin(7rp) vanishes. Thus, we can write 

2m J 2i7i J 



i2Z 



(19) 



(20) 



where the contour of integration goes around the cut of G in the positive direction. 

Other loops can be easily obtained by the methods of and expressed as contour 
integrals. Consider for instance the following case: 




An 8-like contour on the sphere 



The associated Wilson average Wi^2 satisfies —{di + ^2)^1^2 = W1W2, where we note 
Wi for a simple loop enclosing area i and di for derivation with respect to area i (the total 
area A being kept fixed). We now integrate this equation to obtain : 



{2in) 



jdxdy 



^G{x)+G{y)-Aix-A2y 
X + y 



(21) 



Indeed, we have just seen that both sides vary by the same quantity when the areas 1 and 
2 are changed by the same amount; so, we only have to check that they coincide when 
one of the areas vanishes. Suppose A2 = 0: the left-hand side of (^Tj) reduces to Wi, and 
in the right-hand side, we use G{y) ~ ^ for y — > 00 to compute the residue at infinity of 
the y integral, and find Wi again. 

So,(^) holds provided the smallest area is associated with the outer contour. But 
contours can be interchanged: the difference between the integral with y on the outer 
contour and the one with y on the inner contour is given by the (sum of) residue(s) at 
—X (for each value of x), that is 



2i7r J 



[22) 



5 



that is 0, because G is an odd function. We conclude that contours can be interchanged: 
their relative position does not matter, provided they do not intersect (more properly: 
X + y shall not vanish; otherwise, the double integral is still convergent, but it does not 
give the right answer). 

As another example, we consider the Pochhammer contour 




or 




for which the loop equations give: 



(^3 










(-^4 




-di)Wi^2,3A = 




h3W3+4 


(^2 - ^4 


+ ds 


- di)Wi^2,3A = 




^.2+3^^2+3+4 



SO that dsWi^2,3A = ^^1+2+3^^2+3+4 - ^^1+3^^3+4 and 



2,3,4 



(2i7r) 



j) dx dy QG{x)+G{y)-Aix-A4y-A:i{x+y) 



I _ Q-M{x+y) 

X + y 



(23) 



2,4 



(24) 



We take ^3 = to compute -^1,2,4: in this case the Pochhammer contour reduces to an 
8- like loop (^^1,4), and 



(2z7r)2 J 



-A2(x+y) 

X + y 



(25) 



Finally, we have: 
1 



2,3,4 



{2i7T) 



j> dx dy 



J3{x)+G{y)-Aix-A4.y 

X + y 



(1 - (1 - e-^^("+^))(l - e-^^("+^))) (26) 



(here again, the relative position of contours is arbitrary: we only require x + y not to 
vanish) or ^"1,2,3,4 = ^"3+1,3+4 + ^"2+1,2+4 - ^"2+3+1,2+3+4- 



3 The small area phase 

When the total area A = Ai + A2 is less than vr^, the density p obeys a semi-circle law 



A 4 



27r V A 



(27) 
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so that ( ]T8| ) gives: 




This resuh had aheady been obtained in [|I| by a direct estimation of the heat kernel 
0, with windings neglected: 

w{A,,A2) ~ / de,...Nj{{e,-e,fe-'^^^^^^^^^'^ ^Ecos(0O (29) 

At large A^, the angles 9 are distributed according to: 
and finally: 

W{A^,A2) = {cos{d)) = {28) (31) 

Apparently, windings (that is, periodicity for 6''s) can be neglected consistently pro- 
vided the support of p doesn't extend out of [— 7r,7r], which is the long as 
^ + ^ > in particular, one of the areas can be arbitrarily large. However, this 
is not the case with the derivation here given: the total area has to be less than vr^ for the 
semi-circle law to be valid; and, for one large area, Boulatov @ checked the exponential 
law: indeed, he could reproduce the large-area expansion given by Gross and Taylor's 
prescriptions 0. 

A deeper understanding of the contribution of windings seems needed to directly handle 
the heat kernel in these computations. 
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